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Abstract. Observations of up to ten carbon monoxide (CO and isotopomers) transitions are presented to study 
the interstellar C ls O/C 17 ratio towards 21 positions in the nearby (d^-T40pc) low-mass star forming cloud 
p Oph. A map of the C ls O J=l-0 distribution of parts of the cloud is also shown. An average 12 C 18 0/ 12 C 17 
isotopomeric ratio of 4.11 ± 0.14, reflecting the ls O/ 17 isotope ratio, is derived from Large Velocity Gradient 
(LVG) calculations. From LTE column densities we derive a ratio of 4.17±0.26. These calculations also show that 
the kinetic temperature decreases from about 30 K in the cloud envelope to about 10 K in the cloud cores. This 
decrease is accompanied by an increase of the average molecular hydrogen density from 10 4 cm -3 to JjlO 5 cm -3 . 
Towards some lines of sight C 18 optical depths reach values of order unity. 
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1. Introduction 

Abundance ratios of interstellar isotopomers are a power- 
ful tool to study the chemical evolution of the Galaxy. 
One such ratio is that of the rare species of oxygen, 
18 O and 17 O, as measured from the isotopomers of 
CO. For the galactic disk and -center region, Penzias 
(1981) reported average 18 0/ 17 ratios of 3.65±0.15 and 
3.5±0.2, respectively. He also found that the ls O/ 17 ra- 
tio, determined from the integrated line intensity ratios 
/T[ 12 C 18 O(l-0)]di>//T[ 12 C 17 O(l-0)]dw, shows no sig- 
nificant gradient with galactocentric distance R out to 
10 kpc (solar circle: i?o=8-5kpc): the 18 0/ 17 ratios of 
the galactic disk and -center are, within the limits of obser- 
vational accuracy, identical. Models of the chemical evolu- 
tion of the Galaxy by Prantzos et al. (1996) suggest that 
after a few Gyr the ratios in the Galaxy should be in- 
dependent of galactocentric radius. There is, however, a 
discrepancy between the interstellar medium (ISM) val- 
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ues and the much higher (5.5; Anders & Grevesse 1989) 
solar system one. Heikkila, et al. (1998) obtain a low value 
of 1.6±0.3 in the LMC, while in the nuclear starbursts 
NGC253 and NGC494 5 ls O/ 17 - 6.5 (Harrison et al. 
1999; Wang et al. 2004). These results suggest that the 
18 0/ 17 ratio depends on metallicity. 

The sources observed by Penzias are located in a lim- 
ited range of galactocentric radius (and therefore metal 
abundance), and we (Wouterloot et al., in preparation) 
have reobserved these sources with higher angular res- 
olution and have extended our study to sources out to 
i?=16kpc. While Penzias (1981) only observed the J=l-0 
transition, our observations also include the J=2— 1 rota- 
tional lines. The goal of the present paper is to study in 
detail excitation and opacity effects that could affect the 
measured 18 0/ 17 ratios and radial gradients on small 
scales. These effects are usually either ignored or physi- 
cal parameters are derived by assuming a fixed 18 0/ 17 
ratio, so that a more careful study is desirable. 

We have chosen the p Oph cloud (d~140pc) because of 
the large range in column densities found therein, and be- 
cause it is in the solar neighbourhood so that a high linear 
resolution can be attained towards this object. Twenty one 
positions were selected for observations from a C 18 O(l-0) 
map of Wilking & Lada (1983) to have a range of C 18 0(1- 
0) intensities. Towards the positions with the strongest 
lines we not only observed the J=l-0 and 2-1 lines of four 
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isotopomers ( 12 C 16 0, hereafter 12 CO; 13 C 16 0, hereafter 
13 CO; 12 C 18 0, hereafter C 18 0; 12 C 17 0, hereafter C 17 0), 
but we also measured the J=3-2 lines of C 18 and C 17 0. 

2. Observations 

2.1. SEST observations 

Between September 19 and October 5, 1987 we used the 
15-m Swedish-ESO Submillimeter Telescope (SEST) to 
map C 18 O(l-0) towards a part of the p Oph cloud. We 
employed a Schottky receiver in combination with an 
acousto-optical spectrometer (AOS) which had a channel 
separation of about 43 kHz (0.12 kms -1 ). All observations 
were made using frequency switching. The spectra were 
folded and subsequently resampled to a channel width of 
0.24 kms -1 . The rms of the resampled data is typically 

0. 14 K (XX). A 15' x 15' region was observed with a 40" (in 
part of the map 20") raster (the beam size of the SEST at 
110 GHz is about 47"). The mapped region contains the 
cloud cores pOph Bl, pOph C, pOph E, and pOph F, as 
defined from DCO + maps by Lorcn ct al. (1990). 

Between January 26 and 31, 1999 we used two SIS re- 
ceivers at the SEST to observe 12 CO, 13 CO, C 18 0, and 
C 17 J=l-0 and 2-1 towards 21 positions in the p Oph 
cloud. The observed transitions and frequencies are given 
in Cols. 1 and 2 of Table 1. Col. 3 gives the beam size and 
Cols. 4 and 5 the main beam- and moon efficiencies of the 
telescope used (Col. 6). Most intensities in this paper are 
on the T£ scale because then uncertain corrections to Tj£ 
or T m b do not affect the discussed ratios. In some places, 

1. e. when comparing lines from different rotational transi- 
tions, the use of a different scale (main beam brightness 
temperature, T m b, or the average of T£ and T m b) was un- 
avoidable. This is then mentioned explicitly. The pointing 
accuracy was about 5". 

The observed positions were selected from the C 18 
map of Wilking & Lada (1983), and span a large range 
in C 18 intensity and hence in A v ; they are listed in 
Table 2, in order of decreasing intensity (as estimated 
from the Wilking & Lada map). Given are a reference 
number in Col. 1; the position in equatorial coordinates 
in Cols. 2 and 3; the offset positions with respect to 
a(1950)=16 h 24 m 10 s , <5(1950)=-24°23' (this is the aver- 
age of the H 2 CO (Martin-Pintado et al. 1983) and NH 3 
(Zcng ct al. 1984) positions determined for core Bl with 
the 100-m telescope at Effelsberg)[a(2000)=16 h 27 m 11.6 s , 
c5(2000) = -24°29'42"], and the rms in the observed tran- 
sitions in Cols. 6 to 15. 

Towards all positions (1 to 21) we observed each of 
the four isotopomers simultaneously in the J=l-0 and 2-1 
transitions. The observations were made using frequency 
switching and we used the high resolution spectrometer 
(AOS channel-spacing about 43 kHz) split into two equal 
parts. Integration times were chosen to obtain similar sig- 
nal to noise ratios in the C 17 and C 18 spectra in order 
to accurately derive the line ratios. In addition to the sin- 
gle positions we made small (3x3) maps centered on the 
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Fig. 1. Distribution of C 18 O(l-0) emission towards the p Oph 
cloud, a. Integrated between —1 and +8 kms -1 . Contour lev- 
els are 2 to 10 Kkms -1 in steps of 1 Kkms -1 . Small crosses 
indicate the observed positions. Offsets are with respect to 
a(1950) = 16 h 24 m 10 s , (5(1950) = -24°23'. The filled circles are 
those of the 21 positions (see Table 2) that were also observed 
in other isotopomers and that are located within or very near 
the mapped region, b. Peak TJ in the same velocity interval. 
Contour levels are 0.5 to 4.5 K in steps of 0.5 K. The open cir- 
cle indicates the angular resolution of the SEST. The crosses 
indicate in order of decreasing declination the DCO + cores Bl, 
C, E, and F (Loren et al. 1990). 

selected 21 positions on a 20" raster in C 18 (J=2-l) to 
see whether the positions are located in regions with large 
gradients where pointing errors can influence the observed 
line ratios, and to be able to convolve the J=2-l data to 
the J=l-0 angular resolution. 
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SEST 
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-2) 


329330.545 


14" 


0.6 


0.9 


JCMT 


1 "l 


12 C 17 0(3- 
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337061.130 


14" 


0.6 


0.9 


JCMT 





2.2. JCMT observations 

On July 14, 2001, we observed C 17 0(3-2) towards eight 
of the 21 positions with the 15-m James Clerk Maxwell 
Telescope (JCMT) using frequency switching. The ve- 
locity resolution of the autocorrelation spectrometer was 
0.14 kms" 1 and the rms noise level ranged from 0.07 to 
0.14 K, depending on the line intensity (we tried to reach 
similar signal-to- noise ratios at all positions). 

On February 28, 2002 we observed in the same way 
C 18 0(3-2) towards six of the positions previously ob- 
served in C 17 0(3-2). The rms listed in Cols. 8 and 11 
of Table 2 was 0.14 to 0.33 K. 



3. Results 

The C 18 O(l-0) distribution in the mapped region is shown 
in Fig. 1. The map includes the cloud cores pOph Bl, 
pOph C, pOph E, and pOph F (see e.g. the 1.3 mm con- 
tinuum map in Fig. 1 of Motte et al. (1998) which indi- 
cates the locations of the cores A, Bl, B2, C, D, E, and 
F). These cores were originally defined in DCO + maps 
by Loren et al. (1990). We show the emission integrated 
over all velocities (-1 to 8 kms -1 ) in Fig. la and the 
peak TjJ distribution in Fig. lb, respectively. A compar- 
ison between the panels (and Gaussian fits to the lines) 
shows that pOph Bl and pOph C (near offsets (0, 0) and 
(-4,-3) respectively) have relatively narrow lines (about 
1.5 kms -1 ) and a (more) pronounced peak in Fig. lb, 
whereas broader (about 2.3 kms -1 ) lines occur north of 
pOph E, near (—1,-8). At the edges of the map, near 
pOph F, southwest of pOph C, and in between pOph Bl 
and pOph C line widths reach values of about 1.0 kms -1 . 
Compared to the lower angular resolution C 18 O(l-0) map 
in Fig. 2 of Wilking & Lada (1983; beam size 1.1' on a 1' 
or 2' raster), the distribution in Fig. 1 shows finer spatial 
structure. There is also a reasonable correlation between 
the 1.3 mm continuum in Fig. 1 of Motte et al. (1998) 
and the C 18 distribution. The C 18 J=l-0 emission in 
four 1 kms -1 -widc channels is shown in Fig. 2. pOph C 
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Fig. 2. C 18 O(l-0) distribution in four velocity intervals (given 
in the lower right hand corner of each panel). Contour levels 
are 0.5 to 4.5 Kkms - in steps of 0.5 Kkms -1 . 



and pOph E arc mainly observed at 4—5 kms 1 , whereas 
pOph Bl and pOph F show most emission at 3—4 kms -1 . 

All spectra measured towards the first three of 
the twenty-one positions (Table 2) are shown in Fig. 3 
(the spectra towards all positions are published in the 
Appendix in the electronic edition). The velocity inter- 
val is -8 to +16 kms -1 . Most 12 CO and some 13 CO 
spectra show self-absorption (a clear minimum in between 
two peaks in the spectra that is not seen in the lines 
of the rarer CO isotopomers). If flat-topped spectra are 
also considered as an indication for self-absorption, this 
phenomenon occurs at even more positions. The numbers 
in the boxes of the C 18 0(2-l) profiles show the values 
of /TX[C 18 0(2-l)]du for the 20"-spaced nine-point map 
around each position. Some positions show significant (10 
to 15%) intensity gradients where pointing differences be- 
tween C 18 and C 17 (which were not observed simulta- 
neously) could influence the derived line ratios. Equivalent 
line widths (/TXdw/1.06TX[peak]) for C 18 O(l-0) range 
between 0.84 kms -1 (pos. 14) and 1.91 kms -1 (pos. 8). At 
many positions the C 18 and C 17 line profiles show the 
presence of several velocity components, which are more 
pronounced in 13 CO and 12 CO (but at slightly different 
velocities, probably due to self-absorption and saturation 
in the more abundant isotopomers). Two of these veloc- 
ity components are visible (at 3-4 and 4-5 kms -1 ) in the 
C 18 O(l-0) channel maps of Fig. 2. The C 17 spectra also 
show hyperfine structure (see e.g. Lovas & Krupenie 1974). 
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Table 2. Observed sources and rms values (K on a T\ scale) for individual transitions. 



Pos. a(1950) <5(1950) Offset" rms C 17 rms C 18 rms 13 CO rms 12 CO 

hms ° ' " ' 1-0 2-1 3-2 1-0 2-1 3-2 1-0 2-1 1-0 2-1 
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9^ ^4 o 


-24 


97 4^ 


o.oo 


4 7^ 


093 


09Q 


08^ 


0^0 


070 


0.14 


0.12 


90 


90 


93 


2 


16 


23 11.6 


-24 


14 20 


-13.20 


8.27 


0.025 


0.042 


0.140 


0.068 


0.091 


0.21 


0.07 


0.14 


0.27 


0.26 


3 


16 


24 02.1 


-24 


31 33 


-0.53 


-8.55 


0.018 


0.032 


0.071 


0.070 


0.066 


0.19 


0.15 


0.20 


0.30 


0.30 


4 


l fi 


24 02.7 


-24 


29 41 


1 fi7 


fi fiS 


01 Q 


09^ 


081 


0^3 


078 


90 


1 1\ 


0.22 


97 


9fi 


5 


16 


93 98 fi 


-24 


16 34 


-9.43 


6.60 


0.017 


0.031 


0.076 


0.052 


0.077 


0.33 


0.13 


0.20 


0.30 


0.30 


6 


16 


24 13.4 


-24 


33 32 


0.77 


-10.53 


0.026 


0.028 


0.086 


0.060 


0.087 


0.28 


0.17 


0.20 


0.30 


0.33 


7 


l fi 


9^ 31 Q 


-24 


90 nn 


8 fiS 


3 00 


01 fi 


09Q 


1 4^ 


O'iS 


077 




0.14 


0.21 


34 


30 


g 


16 


24 18.1 


-24 


30 21 


1.85 


-7.35 


0.013 


0.020 


0.103 


0.045 


0.051 




0.12 


0.19 


0.35 


0.33 


9 


16 


23 11.2 


-24 


13 30 


-13.38 


9.50 


0.020 


0.029 




0.072 


0.074 
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0.20 


0.34 


0.30 
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1 fi 


£d 4:0. U 


-24 


22 11 




89 


01 Q 


031 




1 04 


0fi3 




0.14 


1 Q 


3Q 


34 


11 


16 


24 26.0 


-24 


27 28 


3.65 


-4.47 


0.023 


0.033 




0.079 


0.059 




0.15 


0.18 


0.33 


0.28 


12 


16 


23 33.8 


-24 


10 48 


-8.25 


12.20 


0.018 


0.025 




0.071 


0.059 




0.14 


0.22 


0.36 


0.33 


13 


16 


23 35.5 


-24 


25 23 


-7.85 


-2.38 


0.012 


0.020 




0.067 


0.060 




0.13 


0.24 


0.38 


0.33 


14 


16 


23 50.0 


-24 


18 20 


-4.55 


4.67 


0.016 


0.026 




0.071 


0.049 




0.16 


0.20 


0.39 


0.36 


15 


16 


23 18.3 


-24 


18 14 


-11.77 


4.77 


0.012 


0.017 




0.076 


0.051 




0.16 


0.23 


0.36 


0.34 


16 


16 


23 20.5 


-24 
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-11.27 


-0.10 


0.0072 


0.014 




0.031 


0.030 




0.15 


0.21 


0.38 


0.35 


17 
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24 20.9 
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15 23 


2.48 


7.62 


0.011 


0.016 
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18 
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09 16 
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0.012 




0.051 


0.041 
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0.32 


0.28 


19 
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-24 
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-14.62 
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0.010 


0.019 




0.038 


0.032 




0.16 


0.21 


0.35 


0.32 


20 
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-24 
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7.62 
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0.034 
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0.13 


0.19 


0.32 


0.38 


21 


16 
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-24 


07 07 


-10.12 


15.88 


0.010 


0.017 




0.042 


0.050 




0.14 


0.23 


0.28 


0.30 



a. With respect to a(1950) = 16 h 24 m 10 s , <5(1950)=-24°23' 



4. Isotopomeric ratios 

4.1. Isotopomeric ratios derived from line intensities 

Integrated line intensities over the velocity interval —1.5 
to 8 kins" 1 (i.e. over all velocity components) and their 
ratios at the measured 21 positions are listed in Table 3 
for C 17 0, C 18 0, and 13 CO. For each position the first 
line gives the intensities (Cols. 2 to 9) and ratios (Cols. 
10 to 14) and the second line the uncertainty therein, ob- 
tained from the rms in the baseline. We note that while 
we obtained the positions from the C 18 O(l-0) map of 
Wilking & Lada (1983) and ordered them according to 
an expected decrease in integrated intensity, we see this 
in our results, but there is significant scatter. This is likely 
caused by a lower angular resolution and undcrsampling 
in the Wilking & Lada map, compared to our data. The 
selected positions cover a large range in C 18 O(l-0) inte- 
grated intensity (from 0.98 to 8.59 Kkms -1 ) (or A v w 
10 to 200 mag (using column densities in Table 4 and 
Frerking et al. 1982)), as intended. 

The resulting C 18 0/C 17 and 13 CO/C 18 integrated 
line intensity ratios for the (1-0) and (2-1) transitions as 
a function of /TJ[C 17 O(l-0)]dw (the most optically thin 
transition) are shown in Figs. 4a-d. Here (not in Table 3, 
following Pcnzias 1981) we corrected the ratios for the 
difference in frequency, which amounts for C 18 0/C 17 
to a factor 1.047 (C 18 0/C 17 0(correctcd) = {u 17 /u 18 ) 2 
C 18 0/C 17 0(observed), see Linkc et al. 1977), but not for 
optical depth and excitation effects. Therefore, if these ef- 



fects are negligible, the corrected ratios should be equiva- 
lent to those of the column densities. Formal errors derived 
from the uncertainties in the line areas (see Table 3) are 
small with respect to the errors introduced by the calibra- 
tion uncertainties. The latter may amount to 7% causing 
an error in the ratios of 10%. The 3x3 C 18 0(2-l) maps 
around the observed positions (see Fig. 3) show that over 
an angular scale of 20" the change in integrated intensity 
is typically 10% or less. This implies that pointing errors 
of 5" cause errors of a few percent in the observed line 
intensity ratios. In Fig. 4b we also show (as filled circles) 
the C 18 0/C 17 ratios derived from the JCMT J=3-2 ob- 
servations. These agree well with the J =2-1 and J =1-0 
results. For comparison we also show in Fig. 4a the average 
result from Penzias (1981) for the galactic disk as a dashed 
line, and the value towards our pos. 1 derived by Bensch 
et al. (2001) from 13 C 18 and 13 C 17 O(l-0) observations 
dotted line. 

At pos. 1, where integrated C 17 line intensities 
and CO column densities are highest, the observed 
C ls O/C 17 ratios are significantly lower than towards the 
other positions (see Figs. 4a, b). This holds for all three ob- 
served rotational transitions. At pos. 10 the C 18 0/C 17 
ratio is lower than at other positions (except pos. 1) 
for J=l but not for ,1=2 1. Omitting pos. 1, the un- 
weighted average of the C 18 0/C 17 integrated intensity 
ratios (including frequency correction) is 3.53±0.48 (sd; 
me 0.11) for the (1-0) and 3.06±0.49 (sd; me 0.11) for 
the (2-1) transition (for average values we derive both 
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Fig. 3. 12 CO, 13 CO, C 18 and C 17 spectra towards the first three of the selected positions given in Table 2. Numbers in 
the C ls O(2-f) panels indicate integrated line intensities in Kkms -1 on a 20" raster centered at the respective position. The 
complete set of spectra can be obtained in FITS format at http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A4-A/ 



the standard deviation (sd) and the error in the mean 
(me=sd/\/]V), which is the most relevant parameter de- 
scribing the uncertainty of the average values). For J=3-2 
(positions 2 to 6) we obtain a ratio of 2.78±0.40 (sd; me 
0.18); for the same positions the J=l-0 and 2-1 ratios 
are 3.37±0.26 (sd; me 0.12) and 2.95±0.18 (sd; me 0.08), 
respectively, suggesting a decrease of the ratio with J, pos- 
sibly because of increasing optical depth (see Sect. 4.2.3). 
At pos. 1 the ratios are 2.18±0.02 (1-0), 1.82±0.02 (2-1), 
and 2.09±0.12 (3-2). 

In Figs. 4c, d 13 CO/C 18 ratios are shown for the (1- 
0) and (2-1) transitions. This ratio may be strongly af- 
fected by 13 C fractionation for small column densities 
(Bally & Langer 1982, Langer et al. 1984), self-shielding 
(van Dishoeck & Black 1988), and by high 13 CO opti- 
cal depths at large column densities, which explains the 
decrease from 10 to 20 in the outer parts of the p Oph 
cloud to about 3 in the cloud centre. This decline in in- 
tegrated line intensity ratios is similar to that seen in 
Barnard 5 (Langer et al. (1989), their Fig. 5), but in 
p Oph 13 CO/C 18 ratios reach even lower values than 
in Barnard 5. Towards the very outer parts of Barnard 5 



where /T[ 13 CO(l-0)]dt; is 1-4 Kkms" 1 , this ratio also 
decreases. In p Oph we do not see this effect, possibly be- 
cause here the 13 CO emission is much stronger towards all 
observed positions. We note that Zielinsky et al. (2000) 
could explain increasing 13 CO/C 18 ratios towards the 
edge of a Photon Dominated Region (PDR) by the pres- 
ence of few big clumps in the cloud center and many small 
clumps at the cloud edge. 

4.2. Isotopomeric ratios derived from LTE column 
densities 

The isotopomeric ratios are likely to be affected by op- 
tical depth effects, which will be strong for CO (most 
positions show self-absorption), significant for 13 CO, and 
not negligible for C 18 0. For C 17 small optical depths 
are expected towards all observed lines of sight. At pos. 
1, which shows the highest C 17 column density, the ob- 
served integrated intensity ratios for J=l-0 C 18 0/ 13 C 18 
and C 17 0/ 13 C 17 are 35.8 and 68.7, respectively (Bensch 
et al. 2001), indicating (assuming that the 12 C/ 13 C ratio 
is about 70; Wilson & Rood 1994) that the C 18 opti- 
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Table 3. Observed C lr O, C ls O, and 13 CO (1-0), (2-1), and (3-2) integrated line intensities and line ratios (i.e. not corrected 
for the frequency difference). For each position the second line gives the formal errors not including calibration uncertainties. 



Pos. 


(1-0) 


C 17 
(2-1) 


(3-2) 


C 18 
(1-0) (2-1) 
Kkms- 1 


(3-2) 


13 

(1-0) 


CO 
(2-1) 


c ls o 
c^o 
(1-0) 


c 18 o 

Cl7o 

(2-1) 


c 18 o 

Cl7o 

(3-2) 


13 co 

C ls O 

(1-0) 


13 co 

C ls O 

(2-1) 


1 


4.12 


3.79 


2.03 


8.59 


6.60 


4.06 


24.45 


15.46 


2.08 


1.74 


2.00 


2.85 


2.34 




0.02 


0.02 


0.10 


0.05 


0.05 


0.16 


0.12 


0.15 


0.02 


0.02 


0.12 


0.02 


0.03 


2 


2.25 


3.52 


3.52 


6.82 


10.07 


8.42 


43.07 


35.38 


3.03 


2.86 


2.39 


6.31 


3.51 




0.03 


0.03 


0.16 


0.07 


0.07 


0.24 


0.08 


0.10 


0.05 


0.03 


0.13 


0.07 


0.03 


3 


2.50 


2.43 


1.22 


8.33 


7.17 


3.13 


27.03 


19.12 


3.33 


2.95 


2.57 


3.24 


2.66 




0.02 


0.02 


0.08 


0.07 


0.05 


0.22 


0.15 


0.15 


0.04 


0.03 


0.25 


0.03 


0.03 


4 


2.62 


2.29 


1.49 


7.52 


6.62 


3.65 


25.33 


18.31 


2.87 


2.89 


2.45 


3.37 


2.77 




0.02 


0.02 


0.09 


0.05 


0.06 


0.23 


0.16 


0.16 


0.03 


0.03 


0.22 


0.03 


0.03 


5 


1.27 


2.69 


3.38 


4.27 


6.74 


8.48 


34.31 


28.68 


3.36 


2.50 


2.51 


8.03 


4.26 




0.02 


0.02 


0.09 


0.05 


0.06 


0.38 


0.14 


0.15 


0.06 


0.03 


0.13 


0.11 


0.04 


6 


2.00 


2.06 


0.95 


7.02 


6.00 


3.21 


22.31 


16.41 


3.51 


2.91 


3.37 


3.18 


2.73 




0.03 


0.02 


0.10 


0.06 


0.06 


0.33 


0.18 


0.15 


0.06 


0.04 


0.49 


0.04 


0.04 


7 


1.97 


3.71 


4.85 


6.31 


8.84 




33.63 


23.46 


3.20 


2.38 




5.33 


2.65 




0.02 


0.02 


0.17 


0.06 


0.06 




0.14 


0.15 


0.04 


0.02 




0.06 


0.02 


8 


1.47 


1.37 


0.57 


5.20 


4.84 




24.02 


15.08 


3.55 


3.53 




4.61 


3.12 




0.01 


0.01 


0.12 


0.05 


0.04 




0.13 


0.14 


0.05 


0.05 




0.05 


0.04 


9 


2.04 


3.44 




7.79 


10.85 




40.61 


29.41 


3.82 


3.15 




5.21 


2.71 




0.02 


0.02 




0.08 


0.05 




0.12 


0.15 


0.05 


0.02 




0.05 


0.02 


10 


1.41 


1.63 




2.97 


4.27 




21.81 


14.73 


2.10 


2.62 




7.35 


3.45 




0.02 


0.02 




0.11 


0.05 




0.14 


0.14 


0.08 


0.05 




0.27 


0.05 


11 


1.25 


1.17 




4.53 


3.80 




20.52 


13.59 


3.62 


3.24 




4.53 


3.58 




0.02 


0.02 




0.08 


0.04 




0.15 


0.13 


0.09 


0.08 




0.09 


0.05 


12 


1.20 


2.36 




4.32 


5.47 




24.53 


16.41 


3.60 


2.32 




5.69 


3.00 




0.02 


0.02 




0.07 


0.04 




0.14 


0.16 


0.08 


0.03 




0.10 


0.04 


13 


0.83 


1.38 




2.40 


3.47 




25.52 


16.69 


2.88 


2.52 




10.63 


4.80 




0.01 


0.01 




0.07 


0.04 




0.13 


0.18 


0.09 


0.04 




0.31 


0.08 


14 


0.95 


1.35 




2.87 


3.30 




16.25 


14.59 


3.00 


2.45 




5.66 


4.43 




0.02 


0.02 




0.07 


0.04 




0.17 


0.15 


0.09 


0.04 




0.16 


0.07 


15 


1.07 


2.06 




3.71 


6.40 




37.64 


33.15 


3.47 


3.11 




10.14 


5.18 




0.01 


0.01 




0.08 


0.04 




0.17 


0.17 


0.08 


0.03 




0.22 


0.04 


16 


0.33 


0.59 




1.11 


1.69 




19.45 


15.78 


3.36 


2.86 




17.55 


9.35 




0.01 


0.01 




0.03 


0.02 




0.16 


0.15 


0.12 


0.06 




0.53 


0.15 


17 


0.65 


0.67 




2.36 


2.56 




16.94 


11.33 


3.65 


3.80 




7.17 


4.43 




0.01 


0.01 




0.06 


0.03 




0.13 


0.16 


0.11 


0.08 




0.19 


0.08 


18 


0.55 


1.05 




1.94 


3.32 




24.13 


17.47 


3.55 


3.15 




12.46 


5.27 




0.01 


0.01 




0.05 


0.03 




0.13 


0.16 


0.12 


0.04 




0.35 


0.07 


19 


0.32 


0.67 




0.98 


1.62 




18.85 


15.80 


3.08 


2.43 




19.20 


9.77 




0.01 


0.01 




0.04 


0.02 




0.16 


0.15 


0.16 


0.06 




0.80 


0.17 


20 


0.37 


0.44 




1.56 


1.78 




17.10 


11.87 


4.20 


4.07 




10.95 


6.65 




0.01 


0.01 




0.03 


0.03 




0.14 


0.14 


0.12 


0.09 




0.26 


0.13 


21 


0.40 


0.88 




1.67 


2.29 




17.92 


12.85 


4.16 


2.61 




10.71 


5.61 




0.01 


0.01 




0.04 


0.04 




0.15 


0.17 


0.15 


0.06 




0.30 


0.12 



cal depth is about 1.5 whereas the C 17 optical depth is 
small. 



4.2.1. Excitation temperature 



Trying to account for optical depth effects, we have 
derived Local Thermodynamical Equilibrium (LTE) col- 
umn densities of 13 CO, C 18 0, and C 17 0. Below we de- 
scribe how we calculate excitation temperatures and op- 
tical depths. 



For the calculation of T cx , both the Tj£ and the T^ b tem- 
peratures scales are relevant. The T m b-scale is strictly 
valid in the extreme case of a source only covering the 
main beam, while the T^-value applies to the opposite ex- 
treme, a very extended source. We have estimated the CO 
excitation temperatures in two different ways, in each of 
which we used both T£ and T^ b . 
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Fig. 4. C 18 0/C 17 and 13 CO/C 18 isotopomeric ratios as a 
function of JTX[C 17 O(l-0)]dw for (a,c) J=l-0, (b,d) J=2-l 
(crosses) and (b) J=3-2 (filled circles). The C 18 0/C 17 and 
13 CO/C ls O line area ratios have been corrected for the fre- 
quency difference (see Sect. 4.1). The dashed line indicates the 
ratio found by Penzias (1981) for the galactic disk, and the 
dotted line denotes the result from Bcnsch et al. (2001) from 
13 C 18 and 13 C 17 O(l-0). 

Firstly, we derive excitation temperatures from the 
peak temperatures of CO(1-0) and CO(2-l) using 



— 5.532 



In 1 



5.532 



(Ti + 0.819) , 
where T w is the peak T[ 12 CO(1-0)] temperature, or 

11.064 



T ex = 11.064 



In 1 



-i -l 



where T21 is 
(method 1). 



the 
This 



peak 
could 



(T 2 i +0.187) 

T[ 12 CO(2-l)] temperature 
underestimate the excita- 
tion temperature whenever there is self-absorption. The 
method helps to constrain excitation temperatures, in par- 
ticular for the outer parts of the cloud. 

Secondly we obtain excitation temperatures from the 
(2-l)/(l-0) ratio of the integrated intensities of C 18 and 
C 17 0, respectively (method 2), assuming that the tran- 
sitions are optically thin (if this is not the case (mainly 
for C 18 0) the excitation temperature will be underesti- 
mated) and that beam filling effects do not affect the ra- 
tio [C 17 0: (2-l)/(l-0) = 4.0 exp(-10.78/T ox ); C 18 0: (2- 
l)/(l-0) = 4.0 exp(-10.54/T ex )]. The J =2-1 C 18 and 
C 17 data were convolved to the J=l-0 resolution using 
the nine point C 18 J=2-l maps described in Sect. 2 (see 
also Fig. 3). The corrections are small (<5% at most po- 
sitions, but 16% at pos. 19). In contrast to method 1 that 
uses optically thick CO lines potentially tracing predomi- 
nantly cloud envelopes, method 2 is based on tracers that 
are more representative of the entire molecular column 
density. A comparison of excitation temperatures derived 
by methods 1 and 2 can show whether there arc temper- 
ature gradients in the cloud. 

Figs. 5a, b show the (2-l)/(l-0) jT^dv ratios for C 18 
and C 17 0, respectively, as a function of fX^[C 17 0(l- 
0)]dv. There is no correlation. In Figs. 5c-f we compare T cx 



derived from J=l-0 (Figs. 5c, e) and J=2 -1 (Figs. 5d,f) 
12 CO T£ (Figs. 5c,d) values (this, method 1, should pro- 
vide upper limits to T ex ) with the excitation temperatures 
of C 18 and C 17 0, derived from the fT^dv (2-l)/(l-0) 
ratios (method 2). We also show in Figs. 5e,f the excitation 
temperatures after converting to T m b using the efficiencies 
in Table 1 (T m b=T£77 moon /?7 rn ,b). The values derived from 
TX[ 12 CO] range from 14.1 K (J =1-0) and 13.0 K (J=2- 
1) at pos. 3 to 32.7 K (1-0) and 30.9 K (2-1) at pos. 
5. These values are for a T£ temperature scale. They are 
even higher for a T^ b scale: the maximum value is 41 K 
at pos. 5. 

Excitation temperatures derived from C 18 and C 17 
(method 2) are generally lower: 6.3 K (pos. 1) to 14.3 K 
(pos. 19) (C 18 0) and 7.2 K (pos. 1) to 21.6 K (pos. 19) 
(C 17 0). Using main beam brightness temperature ratios, 
the values become larger: 7.9 - 26.4 K (C 18 0) and 9.3 - 
66.4 K (C 17 0) for the same positions. 

The OLS bisector mode was used to derive linear re- 
gression coefficients (see Isobe et al. 1990). There is some 
correlation between the excitation temperatures derived 
from 12 CO and C 18 (the correlation coefficient is 0.70 in 
Figs. 5c). For C 17 0, the correlation is slightly less well de- 
fined (the correlation coefficient is 0.66; Fig. 5d). For C 18 
the slope is closest to 1 for the T m b ratios (0.93±0.16), 
while we obtain a slope of 1.97±0.33 for the ratios. 
For C 17 this is reversed: for the T£ ratios the slope is 
1.08±0.15 whereas that for the T mb values it is 0.34±0.08. 
The higher values for 12 CO compared to those from the 
C 18 and C 17 (2-l)/(l-0) ratios can be explained by 
higher kinetic temperatures in the outer parts of the 
clouds (e.g. Castets et al. 1990) from which the 12 CO emis- 
sion mostly originates. 

Because the extent of C 18 and C 17 clumps is most 
likely larger than the main beam, but small compared to 
the size of the Moon we are using in Sects. 4.2.2, 4.2.3, and 
4.3 the average of both efficiencies. The resulting T ex for 
C 18 and C O are compared with each other in Fig. 5g. 
It shows that T cx is well correlated for both isotopomers, 
but for C 17 it can reach higher values than for C 18 0. 
There is no correlation between optical depth and T cx , and 
therefore the lower T ex of C 18 cannot be explained by 
the fact that we did not correct for optically thick J=2— 1 
lines. 



4.2.2. Optical depth 

C 18 optical depths are often derived from C 18 and 
C 17 data by assuming a certain intrinsic C ls O/C 17 
ratio. However the aim of this paper is to determine this 
ratio and therefore this method cannot be applied. 

We tried to fit the C 17 hyperfine components for sev- 
eral positions with small line widths. The optical depth for 
the main hyperfine component of the J=l-0 transition 
had values of less than a few tenths in most cases. The 
fitting is complicated by the presence of more than one 
velocity component, such as a broader underlying com- 
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for pos. 1. This is a little high when considering 18 0/ 17 
ratios of order 4 (Sect. 5) and the results of Bensch et al. 
(2001; see Sect. 4.2). It is equivalent to a value for the main 
hypcrfinc component of about 0.36. At other positions r to t 
ranges from 0.01 to 0.43. For C 18 O(l-0) the optical depth 
is undeterminable (log of negative number) at pos. 1, 3, 4 
and 6 (but see Sect. 4.2). At the other positions it ranges 
from 0.05 to 2.1. 

4.2.3. LTE column density ratios 

Using the excitation temperature and optical depth de- 
j rived above, the column density is calculated from 



5 10 15 20 25 30 20 40 60 

T ex from ratio C 18 0(2-l)/(l-0) T ex from ratio C 17 0(2-l)/(l-0) 



10 - 



from average r\ 



10 20 30 

T ex from ratio C 17 0(2-l)/(l-0) 



Fig. 5. ( 2-l)/(l-0) / T^dv ratios for C ls O (a) and C 17 (b) 
as a function of / TX[C 17 O(l-0)]dv. c,d. T OX [ 12 CO(1-0)] vs T ox 
obtained from the ratio (2-l)/(l-0) for C ls O (c) and C 17 (d) 
using T%_ (crosses) or T m b (squares). e,f. The same as in c,d, 
but for T cx [ 12 CO(2-l)]. In (g) the T cx from C ls O and C 17 
(using the average of Moon and main beam efficiencies) are 
compared. The drawn lines indicate equal values. 



ponent (e.g. at pos. 10, which is also seen in C 18 0), or 
two narrow components (pos. 3, 4, 7). In some cases we 
could fit line widths and velocities to the corresponding 
C 18 O(l-0) spectrum and could use those as input values 
for the fit. Sometimes the fit gave a high optical depth for a 
weak component, which is not realistic. Limited signal-to- 
noise ratios prevented better determinations of the optical 
depth in these components. 

We also derived optical depths from the excitation 
temperatures determined by method 2 using the radia- 
tive transfer equation (see e.g. Rohlfs & Wilson 1996, Eq. 
(14.48)). We find that for excitation temperatures from 
(2-1)/ (1-0) ratios using an average efficiency (as defined 
in Sect. 4.2.1 and Fig. 5g) the highest total C 17 O(l-0) op- 
tical depth (i.e. the resulting peak optical depth if all hy- 
perfine components had the same frequency) is r to t = 1-09 



N, 



_ 3M0 5 1.065 rAv J 



8ir 3 [J, 



J 



-hv/kT m 



with Aw denoting the full width at half maximum intensity 
of the emission line, where the transition is from level J to 
level J — 1; v is the observed frequency, Q is the partition 
function (e.g. Rohlfs & Wilson 1996, Eq. (14.50)), and fi 
the electric dipole moment of the molecule. 

Assuming that the structure of the cloud is somewhere 
in between that of a point source and a very extended 
source, we used here for all following calculations telescope 
efficiencies which are the average of the main beam and 
moon efficiency. Based on the results discussed in the pre- 
vious sections we decided to adopt for C 18 and C 17 
excitation temperatures and optical depths derived us- 
ing the respective (2-l)/(l-0) ratios, and for 13 CO the 
commonly adopted excitation temperature from the cor- 
responding 12 CO transition. In this way we derive two 
sets of LTE column density ratios for the C ls O/C 17 and 
13 CO/C 18 ratios, one based on the J=l-0, the other on 
the 2-1 data. 

The results for the 21 positions are given in Table 4. 
Col. 1 gives the position, Col. 2 the frequency-corrected 
C 18 0/C 17 ratio, and Col. 3 the derived column den- 
sity ratio. The C 17 excitation temperature, total optical 
depth and column density are in Cols. 4 to 6. At each po- 
sition the first row is for J=l-0 and the second row for 
J=2-l. Cols. 7 and 8 give the excitation temperatures for 
C 18 and 12 CO(1-0), respectively. 

The unweighted average N(C 18 0)/N(C ir O) LTE ra- 
tio for the J=l-0 transition (see Fig. 6a) is 4.07±1.32 
(sd; me 0.32), a higher value than that determined by 
Pcnzias (1981), but close to that derived by Bensch et al. 
(2001) from 13 C 18 and 13 C 17 0. Omitting here the high- 
est highly uncertain value (above 8.0, at pos. 11) (the r is 
close to being undetermined and therefore uncertain), the 
ratio becomes 3.81 ± 0.23). At four positions (1,3,4,6) the 
C 18 optical depth was undetermined, but for two of those 
positions (3,6) a value could be derived using the C 17 
excitation temperatures. However we did not use these 
data points to derive the average N(C ls O)/N(C ir O) ra- 
tio. Similarly, for the J =2-1 transition using the same T ex 
as above, the C 18 optical depths were undetermined at 
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Fig. 6. C 18 0/C 17 0((a, b) and la CO/C 18 (c, d) iso- 
topomeric ratios as a function of jTX[C 17 O(l-0)]dt>. The ra- 
tios were derived from LTE column densities using excitation 
temperatures derived from the corresponding (2-1)/ (1-0) ra- 
tios, except for 13 CO where we used T ex derived from 12 CO. 
The filled symbols indicate positions where the C 17 excita- 
tion temperatures were used for the C ls O column density. The 
dashed line indicates the ratio found by Penzias (1981) for the 
galactic disk, and the dotted line marks the result from Bensch 
ct al. (2001) from 13 C 18 and 13 C 17 O(l-0). 



positions 1, 2, 3, 4, 6, 7, 9, and 11, where at positions 
1 to 4 also the C 17 T ex resulted in undetermined r's. 
Ratios derived for this transition are shown in Fig. 6b - 
the average value without the above mentioned positions 
is 4.35±0.35. Omitting here the highest value (above 8.0, 
at pos. 8), the ratio becomes 4.00±0.29. 

V( 13 CO)/iV(C 18 0) ratios are shown in Fig. 6c,d, 
which were derived from the J=l-0 and 2-1 data, respec- 
tively. One can see that after correction for optical depths 
the decrease in ratios towards the cloud center remains. 
This can be explained by real changes in the ratios such 
as fractionation; modelling them is beyond the scope of 
this paper. 

4.3. Isotopomeric ratios derived from Large Velocity 
Gradient (LVG) modelling 

In order to study the excitation of C 18 and C 17 in more 
detail, we also made LVG calculations (e.g. Castor 1970; 
Scoville & Solomon 1974). Simulating the observed line 
intensities, one can estimate column densities, H 2 volume 
densities and kinetic temperatures. The critical point in 
our LTE approach (Sect. 4.2) is the assumption of a sin- 
gle excitation temperature for all transitions of a given 
CO isotopomcr. The LVG calculations provide a way to 
estimate in how far differences in derived column densities 
are accompanied by changes in excitation temperature. In 
principle, such changes could provide line intensity ratios 
that do not directly reflect the 18 0/ 17 isotope ratio and 
an LVG code is a suitable tool to investigate such effects. 

We used the collision rates from Flower (2001) with 
an ortho/para H 2 ratio of 3.0. Taking instead a ratio of 
0.1 (i.e. almost pure para-H 2 ) does not significantly alter 



the results outlined below. Level populations and expected 
line intensities were calculated for 5 K < Tki n < 35 K 
and 10 3 < n(H 2 ) < 10 6 cm~ 3 . For C 18 we assumed an 
average abundance of [C ls O]/[H 2 ] = 1.7 10~ 7 , derived by 
Frerking et al. (1982) for positions in the p Oph cloud. 
In the calculations we used intrinsic ratios of C 18 0/C 17 
between 1.5 and 8.0 in steps of 0.25 or 0.5, and adopted 
a velocity gradient of 5 kms _1 pc _1 , which is appropriate 
for the size and line width of the p Oph clouds. 

Before discussing the results we have to check whether 
the assumptions used by Frerking et al. (1982) to derive 
C 18 column densities, iV(C 18 0), are the same as used 
in the present paper. They observed C 18 O(l-0) towards 
twelve positions and detected eight of them with a max- 
imum Ta of 1.3 K. The positions were selected to have 
relatively low column densities in order to be able to de- 
rive extinctions (A v <10 mag, where depletion (see below) 
is not important). Upper limits to column densities were 
derived assuming that the C 18 excitation temperature 
equals T cx ( 12 CO) and that all levels are populated. To 
derive a lower limit to column densities Frerking et al. 
(1982) assumed the same excitation temperature, but with 
only the J=0 and 1 levels populated. The latter is clearly 
not true because we detected J=2-l emission in excess 
of 1.0 K at positions with [C 18 O(l-0)] < 1.0 K. Then 
Frerking et al. used the average value of upper and lower 
limit. Comparing our LTE column densities (Sect. 4.2.3) 
with values derived using Frerking et al.'s assumptions 
shows that Frerking et al. underestimate column densities 
by a factor of 2. Therefore the abundance is too small by 
that amount. 

There is another reason why the assumption of a con- 
stant abundance is not necessarily correct. Kramer et al. 
(1999) found indications of CO depiction by a factor of 
about 3 in the core of IC 5146 at visual extinction A v of 
about 30 mag. IC 5146 is colder than the p Oph cloud 
and Kramer et al. concluded that C 18 O(l-0) and (2-1) 
are optically thin, which is not the case for all positions in 
p Oph. The depletion starts at visual extinctions of about 
10 mag (Bergin et al. 2002), which would correspond to 
a N(C 18 0) of w2 to 4 10 15 cm" 2 (assuming N(R 2 )/A V 
0.9 10 21 cm _2 mag _1 ; Bohlin ct al. 1978). This is in the 
lower range of the column densities at the observed po- 
sitions (see Table 4). However depletion is not expected 
to affect isotopomeric ratios and we have done some LVG 
calculations using a C 18 abundance which is a factor of 2 
lower than that of Frerking et al. (1982). The results differ 
only marginally from those obtained with an undepleted 
abundance and are therefore inconclusive: at all positions 
the derived H 2 densities are approximately a factor of two 
larger and the kinetic temperatures are slightly smaller. 
Because the effects of depletion and the underestimation 
of column densities by Frerking et al. may balance out, 
we decided to use the Frerking et al. (1982) abundances 
in our calculations. 

Towards the 21 positions the observed and predicted 
values of peak C 18 temperatures (using averages effi- 
ciencies, but also T£ and T m b for some calculations) and 
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Table 4. Derived parameters using average beam efficiencies. 



Pos 


ratio 

C 18 


o/c 17 o 


Tex 


r 

c 17 o 


N 


ex 

c 18 o 


T 

-L CX 

12 CO(1-0) 


ratio 


Tkin 


log[n(H 2 )] 


x 2 




fcorr 


N ratio 
J=\- 
3=2- 


K 

-0 
-1 


(tot) 


cm~ 2 


K 


K 


c 18 o/c 17 o 


K 

LVG 


cm~ 3 




1 


2.18 
1.82 


- 


8.0 


1.09 
1.45 


7.2 10 15 
7.4 10 15 


6.9 


22.9 


4.5 


9.4 


5.19 


4.5 


2 


3.17 
3 no 

o.uu 


4.22 


13.6 


0.19 
n 43 


3.0 10 15 
3 1 i n 15 

O. -L 1U 


12.4 


27.3 


4.0 


17.8 


4.89 


6.8 


3 


3.49 
3 nq 


- 


8.7 


0.30 
n 47 


3.1 10 15 
3 n i n 15 

O.U 1U 


7.7 


15.5 


4.25 


11.5 


4.71 


18 


4 


3.01 
3.03 


- 


7.9 


0.43 
0.51 


3.4 10 15 
3.2 10 15 


7.8 


19.3 


3.75 


11.2 


4.71 


20 


c; 
o 


3 ^9 
2.62 


3 no 

o.uu 

3.76 


24.8 


n 

u.uo 

n i 

U. -LO 


2.2 10 15 
9 3 1 n 15 

^ .O 1U 


1 4 \ 


3(S 3 
ou.o 


3 ^ 
o.o 


23.2 


4 71 


11 


6 


3.67 
3.05 


8.95 


8.7 


0.28 
0.41 


2.4 10 15 
2.4 10 15 


7.4 


18.3 


5.0 


10.3 


4.77 


4.7 


7 


3.35 
2.49 


4.03 
3 ^1 

O.O-L 


16.1 


0.12 
u.oo 


2.7 10 15 
9 q 1 n 15 


10.9 


31.5 


3.88 


15.4 


4.80 


2.7 


8 


3.72 


5.30 


8.0 


0.19 


1.7 10 15 


7.8 


16.6 


4.5 


10.6 


4.47 


2.9 




3.70 


8.56 




0.28 


1.6 10 15 














9 


4.00 
3 3fi 

o.ou 


5.57 
is D3 

u.uo 


13.6 


0.14 
fl 9q 


2.6 10 15 
9 7 i n 15 


10.7 


24.9 


5.0 


14.8 


4.83 


0.9 


10 


2.20 

£j. 1 o 


2.54 
2.91 


10.3 


0.26 
n 38 

u.oo 


1.8 10 15 
i 7 in 15 

-L.I 1U 


12.7 


26.1 


3.0 


12.2 


4.59 


14 


11 


3.78 

3 4fl 


8.11 
q ns 


8.2 


0.23 
0.36 


1.5 10 15 
1.4 10 15 


7.4 


19.2 


5.5 


9.1 


4.71 


0.0 


12 


3.77 
2.43 


3.92 

5.88 


18.4 


0.07 
n 1 7 

U. i- / 


1.7 10 15 
i s in 15 


10.3 


22.0 


4.5 


10.3 


4.71 


14 


13 


3.01 
2.64 


3.12 

3.38 


14.1 


0.08 
0.14 


1.1 10 15 
1.0 10 15 


11.7 


23.0 


3.25 


11.2 


4.41 


0.7 


14 


3.14 
2.56 


4.04 
4.96 


12.1 


0.15 

n 9 c i 


1.2 10 15 
1.2 10 15 


9.6 


23.5 


4.25 


9.7 


4.71 


0.2 


15 


3.63 

3 9fS 


3.51 
3 QS 


18.7 


0.05 
0.12 


1.6 10 15 
i 5 in 15 


15.3 


25.7 


3.88 


18.3 


4.47 


0.5 


16 


3.52 
9 QQ 


3.21 

o.^o 


17.6 


0.03 

u.uo 


4.6 10 14 
4.4 10 14 


13.6 


24.1 


3.5 


11.2 


4.11 


6.1 


17 


3.82 
3 q7 

o.y i 


4.57 
^ 1 fi 


9.3 


0.12 
n 1 q 


7.4 10 14 

71 in 14 


9.5 


19.7 


4.63 


10.8 


4.41 


4.3 


18 


3.72 
3.30 


3.49 
3.77 


17.5 


0.03 
0.08 


7.6 10 14 
7.5 10 14 


14.5 


21.7 


3.75 


18.4 


4.23 


0.5 


19 


3.22 
2.55 


2.36 
2.30 


29.5 


0.01 
0.03 


6.1 10 14 
6.1 10 14 


17.5 


31.3 


3.0 


32.5 


3.87 


9.6 


20 


4.39 
4.26 


4.75 
4.91 


9.5 


0.04 
0.07 


4.0 10 14 
3.8 10 14 


9.0 


19.3 


4.5 


11.2 


3.99 


0.2 


21 


4.36 
2.74 


3.38 
3.74 


23.2 


0.02 
0.07 


6.6 10 14 
6.5 10 14 


11.2 


19.1 


4.25 


11.2 


4.41 


41 



the C 18 0/C 17 ratios for the J=l-0, 2-1 and 3-2 lines 
were compared and x 2 values were calculated using for the 
temperatures an uncertainty a of 10% and for the ratios of 
the three observed rotational transitions a a of 0.15, 0.15, 
and 0.20, respectively. The measured temperatures for the 
J=2-l and 3-2 transitions were convolved to the J=l-0 
beam, as described in Sect. 4.2.1. Cols. 9 to 12 of Table 4 
give the results of the LVG calculations: C 18 0/C 17 ra- 
tio, kinetic temperature, molecular hydrogen density and 



X 2 , indicating the relative goodness of the solutions be- 
tween the different positions. Apparently the scatter in 
isotope ratio is smaller than that obtained with the LTE 
approach. 

The parameters of minimum % 2 are shown in Fig. 7 
(for the assumed C 18 0/C 17 ratio) and Fig. 8 (for T kin 
and log(n(H 2 ))). In Fig. 8b the influence of the assumed 
temperature scale is also shown. The minimum \ 2 values 
range between (perfect fit) and more than 40 (bad fit; for 
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pos. 21), and do not show a systematically lower value for 
either temperature scale. Fig. 8b shows a general decrease 
in density and an increase in T^in with decreasing column 
density (or fXX[C 17 O(l-0)]di;). LVG modelling assumes 
that the density is constant along the line of sight, which is 
probably not true and some lines of sight may have larger 
density gradients than others in the C 17 emitting region. 

For all positions except of pos. 1, Fig. 7 shows a clear 
minimum for some C 18 0/C 17 0. Ratios ranges from 3.0 
to 5.5 with an average value of 4.11±0.65 (sd; me 0.14). 
Giving a higher weight to positions with lower \ 2 increases 
the ratio slightly to 4.21. The exception is pos. 1 which 
has the highest H 2 density, where the \ 2 does not much 
increase for higher ratios. We did not use 13 CO data in 
these calculations because this isotopomer is too much 
affected by fractionation which makes it impossible to as- 
sume a single fractional abundance. Likewise 12 CO data 
show (much) higher excitation temperatures than the de- 
rived kinetic temperatures (see also Fig. 5c-f). This, like 
the increase in Tkin with decreasing column density men- 
tioned above is consistent with kinetic temperatures de- 
creasing towards the cloud interiors from about 20-30 K 
to 10K. 

We have to check what the LVG calculations predict 
for T cx , since LTE assumes that they are equal for all lev- 
els. Fig. 9 shows results of the LVG calculations for an as- 
sumed C 18 0/C 17 ratio of 4.0. In Fig. 9a the dotted lines 
indicate the predicted T CX (J=1— 0) for C 17 0. It is equal 
to T^n at high densities, but for about 5 10 3 < n(H 2 ) < 5 



10 4 



T ex is larger than Tk; n , while below 5 10 cm 



this situation is reversed. The dashed and full-drawn lines 
indicate the difference in T ex ( J=l— 0) between C 17 and 
C 18 0. T cx (C 17 0) is larger than T ox (C 18 0), in agreement 
with the observations (see Fig. 5g). For Tki n < 20 K the 
difference is less than 1 K. At low densities T ex (C ls O) is 
slightly larger than T ex (C 17 0). In Fig. 9b the dotted lines 
show the difference in T cx between the J=2— 1 and J=l— 
transitions of C 17 0. It is small for large n(H2) or low Tki n , 
but can reach values of more than 10 K for Tki n > 20 K and 
n(H2) about 10 cm~ 3 . However the difference of this non- 
LTE effect between C 17 and C 18 (the full-drawn lines 
in Fig. 9b) is much smaller: about 1 K or less in the region 
where most data points are located in Fig. 8. Also these re- 
sults are insensitive to changes in the assumed ortho/para 
H2 ratio: larger changes in the excitation temperatures 
occur at higher temperatures (Tki n >25 K) and lower den- 
sities (n(H 2 )<10 4 cm~ 3 ) than those in the p Oph region. 
This suggests that while LTE calculations will overesti- 
mate column densities both for C 17 and C 18 0, the ratio 
of both column densities will be affected much less. 



5. Discussion 

In calculating the isotope ratios we did not distinguish be- 
tween the different velocity components in the cloud (with 
typical velocity differences of about 2 kms -1 ). The reason 
is the non-Gaussian line shape of many components which 
are confused by the hyperfme structure of C 17 (3 compo- 




ratio C^O/C 1 



Fig. 7. Results of LVG calculations. Plotted is the derived x 2 
value as a function of the assumed C ls O/C 17 ratio for the 
positions 1 to 6. The different positions are distinguished by full 
drawn, dotted, and dashed lines. The filled symbols indicate 
the ratio of minimum x 2 at each of the 21 positions. The circles 
are for the positions 1-6, where J=3-2 data exist. The positions 
associated with the filled squares can be identified in Col. 9 of 
Table 4. The small crosses mark the C 18 0/C 17 ratio used in 
the calculations. 



6 5 



X 

"a 



3 - 



<ti: 



yd ' 



<^3 
'0*^-9 
^3 
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20 



30 



6 - 



6 5 



X 



^ > 



. . . , 



* > 



10 



20 



(K) 



30 



Fig. 8. a. Results of LVG calculations. For each of the 21 po- 
sitions T kin and log(n(H 2 )) are given for the minimum \ 2 dis- 
played in Fig. 7. b. The arrows show how Tkin and log(n(H2)) 
change if T£ (full-drawn arrows) and T^ b (dotted arrows) 
are used. The sizes of the arrowheads are proportional to 
jTX[C 17 O(l-0)]du. In both panels the grid of dots indicate 
the T kin and n(H2) values used for the model calculations. 
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10 



20 



(K) 



30 



Fig. 9. A comparison of the T cx of C 17 and C 18 derived 
from LVG model calculations. The grid of dots indicates the 
T kin and n(H2) values used for these calculations, a. The dot- 
ted lines indicate the non- LTE T cx for C 17 O(l-0) derived from 
the LVG calculations using C 18 O/C 17 O=4.0. Contour levels 
are (left to right) 5 to 50 K in steps of 5 K. The dashed 
and full-drawn lines show the difference in T ex for the J=l— 
transition between C 17 and C ls O: T cx (C 17 0) - T ex (C 18 0). 
Contour levels are -2, -1, -0.5, -0.2, -0.1, -0.05 K (dashed), 

0. 025, 0.25, 0.5, 1, 2, 5, 10 K. Some contour values are la- 
belled, b. The dotted lines show the difference in the non- 
LTE T ox for two transitions of C 17 0: T cx [C 17 O(J=l-0)]- 
T ox [C 17 0(J=2-l)]. Contour levels are (left to right) 0.1, 0.5, 

1, 2, 5, 10, 15, 20, 30 K. The dashed and full-drawn lines 
compare the differences in non- LTE T ex for C 17 and C ls O: 
(T ox [C 17 O(J=l-0)]-T cx [C 17 O(J=2-l)]) - (T ox [C 18 0(J=l)]- 
T ex [C 18 0(J=2)]). Contour levels are -5, -1, -0.5, -0.25, -0.1, 
-0.05, -0.025 (dashed), 0.1, 0.5, 1, 2, 3, 4, 5, 10 K. 



ncnts for J=l-0 and 9 components for J=2-l). Assuming 
a depth along the line of sight of 0.4 pc (10'), the crossing 
time for a velocity of 2 kms -1 would be 2.0 10 5 yr, which is 
probably much less than the lifetime of the cloud. This is 
at least 5.5 10 6 yr, the age of the Upper Scorpius subgroup 
of the Scorpius-Centaurus OB association (de Geus et al. 
1989). Wilking et al. (1989) obtained an upper limit for 
the ages of the T Tauri stars in p Oph of 3.0 10 6 yr. This 
would provide enough time for sufficient mixing of the iso- 
topic constituents of the gas. In addition the types of stars 
which produce 17 O and 18 O are not (yet) present in the 
p Oph cloud (they do exist in the nearby Upper Scorpius 



subgroup of the older Scorpius-Centaurus association), so 
both isotopes are not locally produced and it is unlikely 
that they are enhanced in some parts of the cloud by this 
mechanism (see e.g. Henkel & Mauerberger 1993). It also 
seems unlikely that stellar winds from these associations 
can alter the composition of the cloud significantly. 

The results for the C 18 0/C 17 ratios are summarized 
in Table 5, with the transition and the method used. 

The first three entries in Table 5 indicate that the 
observed ratio is very dependent on the transition used. 
We note that Penzias (1981) corrected his ratios only for 
the difference in frequency and not for optical depth ef- 
fects. This suggests that also the .7=1-0 ratios used in 
Penzias' (1981) galactic study need some more analysis. 
This will be discussed in more detail together with new 
measurements in a forthcoming paper. Within the uncer- 
tainties the average J=l-0 value is equal to the num- 
ber obtained by Penzias (1981) for the galactic plane, 
3.65±0.15. The weighted average ratio of our LTE ratios is 
4.17±0.26, whereas the LVG calculations resulted in a ra- 
tio of 4.11±0.14. Because this method combines all obser- 
vational data we consider this the best result. Bensch et al. 
(2001) detected the almost certainly optically thin transi- 
tions 13 C 18 and 13 C 17 O(l-0) towards pos. 1. Their ob- 
served 13 C 18 0/ 13 C 17 ratio, corrected for the frequency 
difference is 4.23±0.53. Using escape probability models 
they derived a ratio 18 0/ 17 of 4.15±0.52, which is iden- 
tical to (but with a relatively large uncertainty because 
only one position was observed) the presently derived av- 
erage value. Note that the position used by Bensch et al. 
(2001) had to be omitted from our calculations because of 
an undetermined optical depth. 

Recently Ladd (2004) observed C 18 and C 17 O(l-0) 
towards some 600 positions in the Taurus clouds and de- 
rives a ratio of 4.0±0.5, in agreement with the present 
result. However, the C ls O/C 17 ratio appears to de- 
crease with increasing integrated C 17 0(1 — 0) intensity. 
Ladd then concludes that the ratio in the inner parts 
is 2.8±0.4, due to larger self shielding of C 18 in the 
outer parts. The range of integrated C 17 O(l-0) inten- 
sity in Taurus (0.2 to 0.6 Kkms -1 ) is much smaller than 
in p Oph because of the smaller line widths. We think 
that this ratio of 2.8 is not real: our LVG models (us- 
ing a smaller velocity gradient of 2.0 kms _1 pc _1 and 
an intrinsic C 18 0/C 17 ratio of 4.0) can reproduce the 
decrease in ratio with 7(C 17 0) if positions with lower 
J(C 17 0) have a lower density and maybe kinetic temper- 
ature than the points with higher J(C 17 0), which is quite 
possible. In addition, it appears that when self-shielding 
is significant for C 18 0, at N(K 2 ) > 1.22 10 21 (Frerking 
et al. 1982), or A v > 1.3 mag (Bohlin et al. 1978), this 
corresponds to a far-UV extinction by dust of >13 mag 
(Aanncstad & Purcell 1973). This implies that beyond this 
extinction (the (C 17 0) self-shielding H 2 column density is 
even higher) there is too little UV radiation left to affect 
the C 18 0/C 17 ratio, in contradiction with the sugges- 
tion by Ladd (2004). In p Oph a systematic decrease of 
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Table 5. Derived C ls O/C 17 ratios compared with some pre- 
vious results. 



Ratio trans. method omit pos. 
C 18 Q/C 17 Q 

3.53±0.11 J=l-0 freq. corr. 1 
3.06±0.11 J=2-l frcq. corr. 1 
2.78±0.18 J=3-2 freq. corr. 1 
4.07±0.32 J=l-0 from N Ue 1,3,4,6 
4.35±0.44 J=2-l from AT lte 1,2,3,4,6, 

7,9,11 



4.17±0.26 average of J=l-0 and 2-1 
from Mte 

4.11±0.14 J=l-0,2-l,3-2; from LVG 

3.65±0.15 J=l-0 Pcnzias (1981) gal. disk 
freq. corr. 

4.15±0.52 J=l-0 Bensch et al. (2001) 
from 13 C ls O/ 13 C 17 
from cxc. model 

4.23±0.53 J=l-0 Bensch et al. (2001) 
from 13 C ls O/ 13 C 17 
freq. corr. 



C 18 0/C 17 with /(C 17 0) is not seen (see Fig. 4), except 
towards pos. 1. 

6. Summary 

From observations of up to three transitions of C 18 and 
C 17 towards 21 positions in the p Oph cloud we derive 
from LTE and LVG calculations C 18 0/C 17 abundance 
ratios of of 4.17 ± 0.26 and 4.11±0.1, respectively. These 
are expected to be identical to the 18 0/ 17 isotope ra- 
tio. The average molecular hydrogen density towards the 
observed position increases from about 10 cm -3 towards 
the postions with low column densities to 10 5 cm~ 3 to- 
wards positions in the cloud cores. The kinetic tempera- 
tures decrease from 30 K or more at the edge of the cloud 
(as derived from the excitation temperatures of 12 CO) to 
20 K at positions with weak C 17 emission to 10 K in the 
cloud cores. 
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